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Abstract 

A solution is proposed for finding the vector potential and mag- 
netic field of any distribution of currents with axial symmetry. In this 
53 ■ approach, the magnetic field and the vector potential are looked for 

not by solving a differential equation but rather through straightfor- 
ward calculation of integrals of one scalar function. The solution is 
expressed in terms of the associated Legendre polynomials P\ m with 
the index m of the Legendre polynomials assuming one value only, 
m = 1. The solution has the form of a series, with the coefficients of 
the polynomials being combinations of multipole moments. 

Key words: electrodynamics, vector potential, spherical harmon- 
ics, Legendre polynomials, magnetic field. 

1 Introduction 

It is common knowledge that the scalar potential $(r) of any charge distri- 
bution p(r') can be calculated (using the Gaussian absolute system of units) 
by the relation (see, e.g., [lj, Chapter 2) 



In the spherical coordinate system r, (p the <&(r) potential can be con- 
veniently found in many cases by expansion in terms of spherical harmonics 
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Y\ m (see, e.g., [2]). The potential at a point is a sum of actions of all the 
charges involved. Therefore, if we are interested in the potential at a point 
inside a charge distribution, we will have to split the region of integration 
in Eq. (p]) into two parts by a sphere of radius r centered at the pole O. 
As a result, we come to equalities of two types, for r > r' and for r < r' . 
Here r refers to the point with the potential we are interested in, and r' is 
the point over which integration will be performed. Accordingly, the part of 
the potential formed by charges with coordinates r' < r we shall denote with 
$q(i"), and the other part, produced by charges with coordinates r' > r, with 
$G( r )- The expressions for <&q(t) and $o( r ) written in the form appropriate 
for us here can be found, for instance, in [2J and are reproduced below: 



«.M=£t fr>0, (2) 

Z=0 m=—l ' ~r 

where Qi m (r) is a multipole moment of order /,m: 



<Wr) = J jT / / p(r>" Y^', ^dV, (3) 



oo Z 
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<Mr) = EE J^G lm (r)y im (A^ (r < r'), (4) 



Z=0 m= 

where Gi m (r): 



' z V 2Z + 1 



^r) = ^j;jJ^Y^WW'. (5) 

The functions Qi m (r) and Gi m (r) depend on r on the upper and lower 
limits of integration as on a parameter. 

The potential at a given point r is a sum of the potentials calculated for 
r > r' and r < r'\ $(r) = $g(r) + $G(r). 

The spherical functions Yi m can be written as (see, e.g., [2]): 



Y lm (0, <p) = S n 



. (^-H)! , - <m,p f6 x 



In this expression, the integers I, m satisfy the conditions I > 0, — I < 
m < h $m — (— l) m for m > 0, S m = 1 for m < 0. P/ m denotes here the 
associated Legendre polynomial. 
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The vector potential A(r) of any current distribution j(r') in vacuum can 
be calculated from the relation (see, e.g., [1], Chapter 2) 



where c is the electrodynamic constant. 

Relation (j?]) resembles in its form equality (CQ), the difference being that 
integration is performed here not over the scalar but rather over the vector 
function j(r'). If one wanted to employ the above method of potential cal- 
culation with the use of spherical harmonics, one would have to apply all 
the calculations to each component separately, with subsequent vector addi- 
tion of the results thus obtained. This approach may meet with considerable 
difficulties. 

On the other hand, if the currents have axial symmetry, one can suggest a 
fairly simple and straightforward method of calculation of the vector potential 
which reduces essentially to calculating an integral of one scalar function only. 

2 Method of calculation of the vector potential 
of axially symmetric currents 

Assume a system of currents with a volume density j(r') which have axial 
symmetry. Denote the axis of symmetry of the system with Z . Introduce a 
spherical coordinate system r, tp such that the $ = axis coincides with the 
symmetry axis Z. Denote the unit vectors along the r, i?, <p axes by n r , n#, n^,. 
The •§ = 7r/2 plane will be called the equatorial plane. 

The system being symmetric, the currents do not depend on the ip coor- 
dinate and have only one component along it: j(r') = j(r', 'd')n lp i. This means 
that all j(r') are circular currents, and that the vector of any element of the 
current lies in the plane parallel to the equatorial plane (we shall call it the S' 
plane). The symmetry of the current system suggests that the vector poten- 
tial A(r) is also axially symmetric, i.e., the magnitude of the vector potential 
A(r, $) does not depend on coordinate (p. But this means that the magni- 
tude of the vector potential A(r, $) can be determined on any plane passing 
through the symmetry axis Z. The simplest way appears to determine the 
potential on the plane passing through the Z axis and the coordinate p = 0, 
and it is this what we are going to do. Call this plane the M plane. 




(?) 
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Figure 1: 



Calculate the vector potential at a point r lying on the M plane. The 
coordinates of this point are (r, 0). It is convenient for the purpose of this 
analysis to pass through this point a plane S parallel to the equatorial plane. 

Consider our situation in more detail. Vector potential at a point r (lying 
on plane M) is produced by currents flowing throughout the space of interest 
here. Consider the potentials created by the current elements dj(r'). Figure 
1 illustrates our case. 

In Fig. 1, the origin of the spherical system of coordinates is at point O, 
the Z axis points directly at us, and the dashed line visualizes one of the 
circular current trajectories. This trajectory lies in the S' plane. The S and 
S' planes are parallel but do not coincide. The direction of the (p axis at 
point (r, 0) is specified in the figure. We are looking for the potential at 
this point. The (p axis at point (r, 0) is perpendicular to the M plane. The 
r axis lies in the M plane. 

The dashed line L on the 5" plane is perpendicular to plane M, i.e., 
parallel to the (p axis at point (r, 0). 

Figure 1 shows two elements of the circular current with densities rfji(r'i) 
and ^2(^2)- The current element (iji(r'i) produces at point (r, 0) on plane 
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M the potential dAi(r, 0). As seen from vector equality (171) , this potential 
is directed in space in the same way as the current element dji(r\). This 
means that the angles between a direction in space (for instance, straight 
lines perpendicular to plane M) and vectors dji(r'i) and c?Ai(r) should be 
equal (in Fig. 1, these are the p' angles). 

The vector of the current element G?ji(r'i) is confined to plane 5", which 
is parallel to the equatorial plane. As seen from expression (J7J), this means 
that the element of the vector potential rfAi(r) generated by the current 
element rfji(r'i) should likewise lie fully in the plane parallel to the equatorial 
plane, i.e., in plane S (the S and S' planes do not coincide). Therefore, the 
cfAi(r) potential can be resolved into two components, dAi t ±_ and dAx jlp in 
the S plane. These components are shown in Fig. 1: dA\ t ± = dA\ siny/ and 
cL4i i¥ , = dA\ cos ip', where the angle (p' is the coordinate of the current element 



Since, however, the current is axially symmetric, there will always be an 
element of current c?j 2 equal in magnitude to the element dji and located at 
the same distance from point (r, 0) as the dji element. The coordinates of 
this current element are (r', i?', — (p'). It produces at the same point (r, i?, 0) 
the potential dA2(r, 0). Resolving potential dA2(r, 0) in the S plane into 
two components, dA 2> ± and dA-i^, we see immediately that the components 
dAi 7 ± and dA 2t ± cancel, while dAi tV and dA 2ip add (see Fig. 1). Thus, after 
integration over all elements of current dj(r') the potential A(r) will have 
only one component, the one along the ip axis, left. But this implies that 
in determination of the A(r) potential one may, rather than taking into 
account the whole potential dA(r) produced by the current element dj(r'), 
restrict oneself to the projection dA v of this potential onto the p axis: dA v = 
dA cos (p', where <p' is the coordinate of the element of current dj(r'). Said 
otherwise, the potential dA v is generated not by the current element dj(r') 
as a whole but by the projection of this element on the line L parallel to 
vector (p at point (r, 0). 

We can now, on replacing vector j(r') in Eq. ((7j) with a scalar J(r') = 
J (r' , , ip') = j(r' ,$')costp' , calculate the magnitude of the vector potential 



dji(r'i). 




(8) 



The potential at any point r can now be written as 



A(r) = A(r,$,0)n, 



(9) 
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Because expression contains scalar quantities, in calculation of the 
integral one may resort to expanding the integrand in spherical harmonics. 
Note, however, that because expression (jSJ) contains cost^' as a factor, one 
may conveniently use in expressions for Y\ m in place of exponential functions 
employed customarily, trigonometric ones (see. e.g., [3], Chapter 21). In this 
case one will have to replace expression with the following relations: 



21 + 1 


(I 


— m)\ 


2tt 


(I 


+ m)\ 


21 + 1 


(I 


— m)\ 



Pimicos "&) cos m<p, 



P/ m (cos •&) sin mip, 



\ 2tx (1 + m) 
(1 = 0, 1,2,. ..;m = 0,l,2, ...,/)• 
This involves a minor change in the form of expansion of the function. 

For r > r'\ 



(10) 
(11) 



Mr, »,o) = - E (2— + E — + E — J • (12) 

where the multipole moments Qf m and Qf m have the form 

QL(r) = [ j j ArVYLW, v')dV, (13) 

QL(r) = J J JVVYLW, tp'W- (14) 

For r < r': 



1 CO / A I -1 t £ \ 

A G (r, 0, 0) = - £ J ^r— rr'G&y^ + £ ^LYL + E rlG im Y im , 

1=0 ' \ Z m=l m=l / 

(15) 

with the multipole moments Gf m and Gf m 
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our) = ^ril^y^wyv, (i6) 

The function cosy/ entering the expression for J(r') also belongs to the 
system of functions Yf m (see expression (fTUj) ). Therefore, because of the func- 
tions cos my?' and sin my/ being orthogonal, part of the multipole moments 
will vanish in the course of integration of expressions ( 1T31) , (TT4l) and ( IT6l) , ( TTT|) 
(orthogonality in index m). The only terms left will be those with m = 1, 
i.e., the terms Qfi{r) and G^(r). As a result, the double series (TT2I) and (TT5j) 
will become single series and look as 

A^,0) = \t^f^ forr>r', (18) 

1 °° / ZL 7T 

A G (r,tf,0) = for r < r'. (19) 

The next step appears to be substitution in these series of the expres- 
sion for Y L \ from formula f llOp . We have to remember, however, that we are 
looking for the potential on the ip = plane only. For all other values of 
the (p coordinate, the relations we have derived will yield wrong values of the 
potential. Therefore, prior to substituting the expression for Yfi from relation 
(fTUj) . we will have to zero the coordinate (p. We obtain naturally cosy = 1. 
Now the potentials take on the form 

M^o)4g^MM« forr> „, (20) 

where the multipole moments are 

Qf 1 {r)= [ r T j(r\$y' l Pi 1 (cos$ l )r' 2 sm'&'d'&'dr', (21) 
Jo Jo 

1 °° 2-7T 

AgMO) = -ETTT^T^nW^Ucos^) for r < r', (22) 
c ^1(1 + 1) 

with the multipole moments G{[ 
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Gf x (r ) = f°° r 3{r '^ P a (cos $')r' 2 sin V'M'dr'. (23) 

Jr Jo r' i+L 

In these formulas, j(r', is the circular current density, and (J/i(r) and 
G{^{r) are the parts of the multipole moments left after integration over the 
ip coordinate (truncated multipole moments). The vector potential at any 
point in space can be written as 

A(r) = [Aq(t, 0, 0) + A G (r, 0, 0)K, (24) 
whence we finally come to 

Mr) = \ t + r'GKr)) P,(cos^)n,. (25) 



This is the final expression for the vector potential A(r). As evident from 
Eq. (I25p . vector potential A(r) has only one component, A v . 

Now we impose an additional constraint, namely, that the currents are 
symmetrical with respect to the equatorial plane (which quite often is the 
case). In this case, expression (125]) may contain only functions Pn symmetric 
about the equator, i.e., P2n+i,i functions. Now expression fl2"5]) for the vector 
potential has to be replaced by 

A « - \ £ > + l)(2 n + l) + ^ G £, +mW ) Wco.*)* 

(26) 

It may be of interest to consider a particular relevant case. Consider 
axially symmetric currents with a sin $ dependence on the $ angle. This case 
becomes realized, for instance, if the currents are generated by rotation of 
charges, with all the charges in a layer of thickness dr distributed uniformly 
over the angle d and rotating with the same angular velocity. 

The sin $ function belongs to the system of associated Legendre functions. 
Therefore, the Pn functions being orthogonal with respect to index /, the 
larger part of the terms in the Qai r ) an d G{[(r) expansions (see expressions 
(12T|) and (123]) ) will vanish in calculation, with only Qfi(r) and G^\(r) terms 
with index / = 1 retained. In this case, the potential will come out as a sum 
of two terms only: 
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A(r) 



* (QUr) 



c \ r z 



+ rG^ 1 (r) sin^n^ 



(27) 



3 Magnetic fields of axially symmetric currents 

Magnetic field H(r) can be calculated with the general expression H(r) = 
V x A(r). The vector potential of axially symmetric currents has only one 
component, A^: A(r) = A(r, i?)]!^. With this in mind, and using the general 
expression for V x A(r) in the spherical coordinate system (see, for instance, 
P], Chapter 6) 



V x A 



1 ( 9 (dn&A ) dA * 
rsint? yd'd v dip 



rsm.'d dip 
we finally come to 



1 dA r 1 / d 



r \dr 



(rA, 



1 d{rAt 



I dA r 



r dr 



r d'd 



n 



tt „ * A COS?? . 1 d . \ (I 

H = V x A = -— ^ + -— A v n r - - 

\ r sin v r dv \r 



K + (28) 



with A<p expressed by formula (125 



Consider the expression jrAu, in any of the terms of expansion (12511 : 



d 



A, 



1 2vr 



(/ + l) Qf x + Jr^Gg + + r^Gg j Pn(cos^). 

(29) 



„Z+2 



J A. 

<9r 



dr" v c 1(1 + 1) \ 
By the theorem of Leibniz-Newton (see, e.g., [3], Chapter 4, item 4.6-5) 

(/ "' f(l)<ll /(,■). 



the last two terms in expression ([29]) cancel, with the final expression taking 
on the form 
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Hfr) 



E 



2tt 



-i(J + l) 



r l+2 



.1-1 



'costf dPn{cos'd) 
-Pixicoav) + 



sin?? 



n r + 



(30) 



r l+2 



-(I + l)r l - 1 Gf 1 (r)\P ll (cos'&) 



If the currents are symmetric about the equator, expression ( 1301) . similar 
to formula (126]) . will contain only equatorially symmetric functions P2n+i,i> 
and summation over n should be performed from to oo. 

Consider a few particular cases. In the case covered by potential ( |27|) . the 
magnetic field should be written as 



H 



7T 

C 



'Qti 



+ G\\ 2 cos n r + 



li 



sin $11,9 



(31) 



To illustrate the use of this expression, consider the following simple prob- 
lem: rotation of a charged sphere about one of the diameters. Let the sphere 
of radius R, with charge q distributed uniformly over its surface, rotate about 
the Z axis with an angular velocity w. Calculate the magnetic field and the 
vector potential inside and outside the sphere. 

The density of the surface current produced by rotation of the sphere can 
be expressed in the following way: 



qw 



5{r — R) sin-^n^. 



(32) 



Substituting the expression for current density ( 1321) in relations ( 1161) and 
(fT7|) . we come to the following expressions for the (truncated) multipole mo- 
ments 



Q n — 0, G u 



qw 
3^R 



Qn 



qwR 2 
3vr 







for r < R, 



for r > R. 



Now the magnetic field and the vector potential acquire their final form 
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H = — ^(cos$n r — sin #11$), A = sin^n^ for r < R, (33) 
3cR 3cR 



2qwR 2 , . 1 . n \ » qwR 2 . „ 

H = — (cosi/n r H — smi/n^J, A = — sinxm^ tor r > R. 

3cH 2 3cr 2 

(34) 

Significantly, this technique permits one to determine the magnetic field 
without prior calculation of the vector potential. 

The solution to this problem can be found, for instance, in monographs 
[lj (problem 2.85) and [2j (problem 253): 



2qw „ 3r(m ■ r) m 

H = for r < R, H = — — — '- - — for r > R, 

3cR r 5 r 6 

where m = ^^rw is the magnetic moment of the system. One can readily 
verify that these expressions are identical. In monographs [T] and [2] it is 
proposed to solve this problem by differential calculus. Our approach, by 
contrast, has yielded the answers to it by straightforward calculation. 
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